Disclination classes, fractional excitations, and the melting of quantum liquid crystals 
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We consider how fractional excitations bound to a dislocation evolve as the dislocation is separated 
into a pair of disclinations. We show that some dislocation-bound excitations (such as Majorana 
modes and half-quantum vortices) are possible only if the elementary dislocation consists of two 
incquivalent disclinations, as is the case for stripes or square lattices but not for triangular lattices. 
The existence of multiple inequivalent disclination classes governs the two-dimensional melting of 
quantum liquid crystals (i.e., nematics and hexatics), determining whether superfluidity and orien- 
tational order can simultaneously vanish at a continuous transition. 

PACS numbers: 



Spatial order (i.e., crystallinity) is found intertwined 
with exotic quantum-mechanical order in various sys- 
tems: for instance, with unconventional superfluidity, 
in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
and related states [IH7]; and with topological electronic 
structure, in weak and crystalline topological insula- 
tors/superconductors [HHH]- In such cases, the de- 
fects of the translational crystalline order — i.e., disloca- 
tions — can bind excitations or defects of the intertwined 
quantum-mechanical order, such as Majorana bound 
states, protected helical edge modes, or half-quantum 
vortices [THE)]. More generally, such composite defects, 
and their proliferation, are central themes in the theory 
of deconfmed quantum criticality |17j . A common fea- 
ture of the dislocation-bound states noted above is that 
they are 11 Z 2 " in character, i.e., two of them combine to 
give a conventional bound state which is trivial in some 
sense. 

The fact that dislocations can bind such Z 2 states 
raises a conceptual puzzle because dislocations are them- 
selves composite defects: a dislocation is a bound state 
of two disclinations, which are topological defects of ori- 
entational crystalline order [T8H22] . When a dislocation 
is split into two widely separated disclinations (Fig. [IJ , 
there is prima facie no straightforward way for the ac- 
companying bound state to divide or to stretch between 
the two disclinations; nor is there an obvious principle 
for assigning the bound state to either disclination. In 
this work we immediately present the resolution to this 
puzzle by explaining why the bound state attaches itself 
to one of the two disclinations, and then we go on to 
discuss the non-trivial consequences of this concept. 

Following the work of Ref. [23], we note that cer- 
tain lattices (including stripes and square lattices) have 
topologically inequivalent classes of disclinations, some 
of which support exotic bound states on their own. In 
fact, a dislocation carrying a Z 2 bound state is always 
composed of one trivial and one nontrivial disclination 
and this composition determines how the bound state is 
resolved when a dislocation splits: it attaches itself to the 



non-trivial disclination. An immediate corollary is that 
in lattices without multiple, inequivalent types of discli- 
nations, there will be no non-trivial Z 2 bound states on 
dislocations either. For the case where translational order 
is intertwined with superfluidity, we show that the pres- 
ence, or absence, of multiple disclination classes modifies 
the structure of the phase diagram in two dimensions. 
For example, in lattices that have only one class of el- 
ementary disclination (e.g., hexagonal lattices [21]), su- 
perfluid and liquid-crystalline order can generically van- 
ish at once through a continuous transition; in lattices 
with inequivalent disclinations, however, this is forbid- 
den except possibly at a multicritical point. While this 
work only explicitly considers superffuid order, the ar- 
guments and main result will also apply to cases where 
dislocations carry other, perhaps magnetic, defects (as in 

Refs. nam). 

The predictions of our analysis are straightforward to 
test in ultra-cold atomic experiments with FFLO liq- 
uid crystals or with spin-orbit coupled condensates [2"T1 - 
133] . Furthermore, our central result regarding inequiv- 
alent disclinations has immediate experimental conse- 
quences for weak topological insulators and superconduc- 
tors where our arguments imply that it is the disclina- 




FIG. 1: (a), (b) ±7r/2-disclination dipoles that fuse into dis- 
locations with odd, even Burgers' vectors respectively. 
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FIG. 2: (a) Two inequivalent classes of disclinations of a two- 
dimensional striped superfluid. Lines denote nodes of the 
order parameter; the sign of the condensate between these is 
indicated on the left. The crest-centered disclination (red cir- 
cle) involves no sign-frustration, whereas the trough-centered 
disclination (red square) does involve sign- frustration, (b), 
(c) Examples of the two disclination classes for a square- 
lattice striped superfluid, indicating the absence and presence 
of sign-frustration respectively. 



tion core-size, rather than the dislocation core-size, that 
sets the confinement scale of bound states, hybridization 
energy gaps, etc. This prediction is particularly salient 
to experiments, as the conditions under which disloca- 
tions are common and easy to probe (i.e., those under 
which crystalline order is weak) are precisely those un- 
der which the dislocations consist of two weakly bound 
disclinations. 

General classification of disclinations. We begin by 
reviewing the results of Ref. [23]; we focus on three 
cases, namely stripes, square lattices, and triangular lat- 
tices. A crucial observation in the analysis of Ref. [23 is 
that the holonomy around a disclination can take mul- 
tiple discrete values depending on where the disclina- 
tion is centered. The holonomy h is simply the element 
in the planar space group accumulated along a loop in 
the lattice around the disclination. It contains a rota- 
tion component (Frank angle) and translation compo- 
nent (Burgers' vector) [TB]. The group element is path 
independent up to conjugation h M- ThT^ 1 by a trans- 
lation T of the starting point of the loop. Disclinations 
are therefore topologically classified by conjugacy classes 
of the planar space group; these classes correspond to 
point-group rotations with inequivalent rotation centers. 
Thus, stripes have two classes of ir disclinations, namely 
disclinations centered on a stripe or halfway between 
two stripes [Fig. |2|a)], and an elementary dislocation 
contains one disclination of each class. For square lat- 
tices, there is a similar Zi classification of 7r/2 disclina- 
tions [Fig. [2jb),(c)], which are either vertex-centered or 



plaquette-centered, and again an elementary dislocation 
contains one disclination of each class. However, for a tri- 
angular lattice there is only one class of 7r/3-disclinations 
and 27r/3-disclinations are .^-classified. One can easily 
check that, because 2 and 3 are coprime, none of the three 
types of 27r/3 disclinations carries a Zi bound state, and 
it follows that all dislocations are topologically trivial. 

Striped superfluids. We now address the nature of 
disclination- vortex bound states, beginning with the case 
of a striped superfluid that is described by a Bose conden- 
sate of the form ?/>(x) = ip cos(k • x + <f>) exp(i#), where ip 
is taken to be a fixed amplitude and 9 and 4> are phases 
associated with the superfluid and translational Gold- 
stone modes respectively. The canonical example of such 
a state is the Fulde-Ferrell-Larkin-Ovchinnikov state pro- 
posed for spin-polarized superconductors |3H [35] ; a 
closely related state, which may occur in layered cuprate 
superconductors, is the pair-density-wave state [5j. Fur- 
thermore, ultracold atomic gases with spin-orbit cou- 
pling offer a purely bosonic analog of these states [2"51 - 
I3UI321I3S]- For specificity, we address the rotationally in- 
variant case; thus, our analysis directly applies to FFLO 
states and to Bose gases subject to a Rashba spin-orbit 
coupling [28 30 , 33 . As we are primarily concerned with 
topological features rather than energetics we could re- 
store weak anisotropy without affecting our conclusions. 

The striped superfluid has three easily identifiable 
types of topological defects [5J: (i) a pure dislocation, in 
which <f> winds by 27r, (ii) a pure vortex, in which 9 winds 
by 27T, and (iii) a half-vortex-dislocation bound state, in 
which both <f> and 9 wind by ir. Since defects of type (i) 
and (ii) can be written as dipoles of (iii), it suffices to 
understand how the vortex-dislocation bound state sep- 
arates into disclinations. The Z 2 classification of discli- 
nations for stripes given above provides a direct (and 
hitherto undiscussed) geometric criterion for understand- 
ing how this separation occurs. We know that the half- 
quantum vortex bound to a dislocation compensates for 
the change in sign of the order parameter going around 
the dislocation. Of the two kinds of possible disclinations, 
only the node-centered disclination (type-A) involves a 
frustration of the order-parameter sign. Therefore, as an 
elementary vortex-dislocation bound state consists of one 
node-centered and one anti-node centered (type B) discli- 
nation, the type-A disclination inherits the half-vortex 
when the dislocation unbinds. 

We now turn to the melting of a striped superfluid. As 
is well known [3], in a rotation-invariant two-dimensional 
system, dislocations proliferate at any nonzero temper- 
ature, reducing the smectic to a nematic. While cold 
atom systems are never exactly rotation symmetric, the 
temperature scale set by the anisotropy (at which disloca- 
tions proliferate) can be taken to be much lower than that 
at which the nematic melts. In what follows, therefore, 
we consider this nematic phase. From above we know a 
type-A disclination carries both a Frank angle (an orien- 
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tational "charge") and a vorticity (a superfluid charge), 
whereas a type-B disclination only has an orientational 
charge. In what follows we denote the elementary discli- 
nations and their composites using the notation (a, j3) 
where a is the orientational charge and (3 is the superfluid 
charge. In this notation, a pure dislocation is (0, 0) (and 
can be constructed as a dipole of two type-A or type-B 
disclinations with opposite a); a pure vortex has charge 
(0, 2) and can be composed of two type-A disclinations of 
opposite a (Frank angle); and finally, a half- vortex half- 
dislocation bound state consists of one type-A and one 
type-B disclination, and thus carries charges (0, 1). 

The low-energy theory of a nematic superfluid is given, 
in the one-constant approximation |18) , by the expression 

H el . = k [(V-n) 2 + (Vxn) 2 ] + Ps \V0\\ (1) 

where n is a unit vector normal to the local orientation 
of the nematic, k is the Frank constant governing orien- 
tational stiffness, and p s is the superfluid stiffness. This 
elastic theory can be analyzed using standard Coulomb 
gas methods [37] ; this involves transforming the elas- 
tic theory to a theory in terms of the topological defects 
discussed above: 

n c = £lW«l 2 + ^|v0 s | 2 (2) 

- X! 9a,fi cos(7ra0„) COS(7T/30 S ). 

a/3 

One can now compute the temperature at which a given 
defect unbinds by computing the scaling dimensions of 
the g a fi of every topological defect (a, j3) around the non- 
interacting fixed point defined by the quadratic terms. 
These scaling dimensions are generically given by the ex- 
pression, A Q/ 3 = {n/T)[a 2 K + /3 2 p s ] [5]; a defect with 
indices (a, /3) proliferates when A Q p = 2. 

In the nematic superfluid, all charge is confined, so that 
the basic components of the equilibrium state are neutral 
dipoles of two type-A or two type-B disclinations. The 
type-A disclinations are bound by both superfluid and 
elastic charge; thus, they are more tightly bound than 
the type-B disclinations, in this phase, and can never 
proliferate before the type-B disclinations. Thus, from 
our analysis we find that there is generically no direct 
transition from the nematic superfluid to an isotropic 
normal phase contrary to Ref. [38] • In the Coulomb gas 
language, one can straightforwardly see this as follows. 
Depending on the values of p s /n, the excitations with 
the lowest threshold for unbinding are type-B disclina- 
tions [(1,0)] or oricntationally neutral half-vortex half- 
dislocations [(0, 1)]; other defects are always less relevant 
than these. For small superfluid densities, the half-vortex 
half-dislocation is the first to proliferate, leading to a 
nematic but non-superfluid state. For larger superfluid 




FIG. 3: (a) Phase diagram of a striped, or square-lattice, su- 
perfluid in two dimensions, for a fixed Frank constant, as a 
function of superfluid stiffness and temperature. Solid lines 
denote phase transitions, due to the proliferation of defects 
(indexed by their orientational and superfluid charges); the 
dashed line represents the temperature at which type-A discli- 
nations would (hypothetically) have proliferated if the sys- 
tem were still ordered. The four phases are nematic/hexatic 
(Ne/He), nematic/hexatic superfluid (Ne/He-SF), isotropic 
superfluid (SF), and isotropic normal (I), (b) Phase diagram 
of a triangular-lattice superfluid [case (a)]. The notation is 
the same as that in (a); note that a direct transition from the 
hexatic superfluid to the isotropic normal state, absent in (a), 
appears here. 

densities, the type-B disclination is the first to prolif- 
erate, giving rise to an isotropic superfluid state. The 
resulting phase diagram is shown in Fig. [3](a). At this 
level of analysis the two transitions meet at a multicrit- 
ical point, at n = p s , which is the only point at which 
a direct transition out of the nematic superfluid into a 
normal phase is possible; however, nonlinearities not in- 
cluded in this analysis might change the phase diagram 
in the vicinity of the multicritical point. (Note that, in 
experiments with spin-orbit coupled Bose gases, one can 
tune the ratio n/p s by tuning the spin-dependent scat- 
tering lengths: n/p s <C 1 near the stripe-to-plane- wave 
phase transition [29 "33], while the opposite limit holds 
deep in the striped phase.) 

We now extend our analysis to the case of hexagonal 
and square-lattice pair-density waves with sign-changing 
order parameters. These can arise due to nonlinear in- 
teractions between stripes, or can be directly engineered 
by tuning interactions (e.g., cavity-mediated [321 E0] or 
fermion-mediated interactions) to be peaked at certain 
momenta. A specific case of a pair-density-wave state 
on a hexagonal lattice was recently discussed in Ref. \7\ 
with an analysis that ignored the effects of disclinations. 
However, disclinations play a crucial role in the melting 
transition and we find, remarkably, that the nature of 
two-dimensional melting depends strongly on the lattice 
geometry, via the classification of disclinations. 

Unlike striped states, crystalline states can sustain 
quasi-long-range order at finite temperatures in two di- 
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mensions; however, they generally melt in two steps, with 
the dislocations proliferating first and then the discli- 
nations [19H21| 141] . We shall analyze the intermediate 
"hexatic" phase, in which pure dislocations have prolif- 
erated (so that translational order is short-range) but 
disclinations are still confined. For square lattices (i.e., 
checkerboard arrangements of the order parameter), as 
with stripes, there are two inequivalent classes of discli- 
nations; as can be seen in Fig. [2j only one of the two 
classes has sign-frustration and thus binds a half-vortex. 
Thus, the overall phase diagram in the square-lattice case 
is exactly as in the striped case: the first defect to pro- 
liferate is either the neutral disclination or the vortex- 
dislocation bound state, and in either case the transition 
out of the hexatic superfluid is into either a hexatic or a 
superfluid state. 

We now turn to triangular/hexagonal lattices where 
we must distinguish between two types of order param- 
eter configurations: (a) that in which the nodes of the 
order parameter form a triangular lattice [Fig. [4ja)] , and 
the antinodes occupy the dual honeycomb lattice; and 
(b) that in which the antinodes form a triangular lattice 
[Fig. |4jb)], and the nodes occupy the dual honeycomb 
lattice. This complication does not appear in the square 
lattice because it is self-dual. We consider case (a) first. 
We see the ±7r/3 disclinations on the triangular lattice in 
Fig. [5^e) break the •, o-sublattice order that corresponds 
to the opposite signs in Fig.[4|a) and therefore invariably 
carries a half-vortex, i.e. it carries charges (1,1). These 
elementary disclinations can be combined to form pure 
vortices (0, 2) or pure double disclinations (2, 0) equiv- 
alent to the 27r/3 disclinations in Fig[5]ja-d). Depend- 
ing on p s /K, any of the three defects might proliferate 
first. In particular, the proliferation of elementary (1, 1) 
disclination-vortex states, which is favored when the ori- 
entational and superfluid stiffness are comparable [i.e., 
1/4 < n/p s < 4], can give rise to a direct Kosterlitz- 
Thoulcss transition between the nematic superfluid and 
an isotropic normal state, as indicated in Fig. [3jb) . (The 
formal structure of the Coulomb-gas problem, and thus 
the overall structure of the phase diagram, are then iden- 
tical to that in Ref. [5].) 

We now briefly comment on case (b), in which the 
antinodes form a hexagonal lattice [7]; one can alterna- 
tively regard this state as a honeycomb lattice in which 
the »-sites are occupied by vortices and the o-sites by 
antivortices in Fig [5j Formally, this situation is some- 
what different from that discussed so far, as the rela- 
tive phase between neighboring antinodes is not n (i.e., 
a minus sign) but is instead 27r/3. This modulation of 
the superfluid phase can be understood as a tri-coloring 
structure in the honeycomb plaquettes in Fig [5] where no 
adjacent plaquettes carry same color. It was argued in 
Refs. [3 [HI H3] that such a superfluid supports a topo- 
logical defect with a phase twist of 27r/3 (i.e., a vorticity 
of 1/3) bound to an elementary dislocation. In the hex- 
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FIG. 4: Two possible varieties of hexagonal-lattice superfluid. 
(a) Hexagonal lattice of nodes of the order parameter, corre- 
sponding to a honeycomb lattice of the particle density, +, — 
represent the sign of Bose condensate ip; (b) Hexagonal lattice 
in the particle density, corresponding to a honeycomb lattice 
of vortices, 0, ±27r/3 are phases of the condensate ip. 
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FIG. 5: (a-d) Inequivalent trivial ±27r/3-disclinations. None 
breaks the a, o-sublattice order (no nearest vertices are of 
same type). Tri-coloring ordering (no adjacent plaquettes 
have same color) is preserved for (a) and (b) and is broken for 
(c) and (d). (e) Topological ±7r/3-disclination dipole. Each 
breaks •, o-sublattice and tri-coloring order. 



atic phase, dislocations are screened; hence this defect 
carries only the elementary superfluid charge (i.e., it is 
simply an elementary vortex). When the orientational 
stiffness exceeds the superfluid stiffness, this proliferates; 
in the opposite limit, the first defects to proliferate are 
the plaquette-centered 27r/3 disclinations [Fig. J5ja,b)], 
which do not violate the tri-coloring pattern, and are 
therefore neutral with respect to superfluidity. Thus the 
overall phase diagram has the form shown in Fig. [3ja) . 
Note that such a phase diagram could arise on the trian- 
gular lattice precisely because the superfluid defect was 
Z 3 , rather than Z 2 , in character. A nontrivial Z 3 defect 
can bind to a dislocation even if the dislocation consists 
of two equivalent disclinations, because 2 = 14-1 and 
1 = 2 + 2 mod 3, so that the original puzzle discussed 
in the introduction does not apply. Thus, case (b) of the 
triangular-lattice superfluid can be regarded as a simple 
counterexample illustrating the crucial role played in our 
analysis by the Z2 nature of the exotic bound state. 
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Finally, we sketch the implications of our results for 
weak topological insulators and superconductors. It is 
well-known that the presence of low-energy bound states 
on a dislocation is determined by the weak topological in- 
dex |15j ; as the weak-index is always zero for a triangular 
lattice, it follows that there are no bound states. On the 
square lattice, one can explicitly see from the construc- 
tion of Ref. [23] that a weak topological invariant entails 
the presence of inequivalent disclinations. We note that, 
although no analogous classification of topological insu- 
lators (i.e., not assuming particle-hole symmetry) has yet 
been performed, one can explicitly see that in the case 
of a weak topological insulator layered along the z axis, 
an edge dislocation line along the z axis containing he- 
lical edge modes can be decomposed into a disclination 
line that carries helical modes and one that does not. 
These observations imply that the transverse size of the 
bound states is set by that of the disclination core, which 
is typically smaller than the dislocation core. The pres- 
ence of bound states also generates an effective attraction 
between disclinations of the same topological type; how- 
ever, one expects this to fall off exponentially with sepa- 
ration and therefore not to have qualitative implications 
for melting transitions. 

In summary, we have shown that the presence, or ab- 
sence, of inequivalent disclination classes has a qualita- 
tive impact on the phase diagrams of pair-density-wave 
states: if there arc inequivalent disclination classes, the 
nematic superfluid must melt in two steps; whereas if all 
disclinations arc equivalent, then a direct transition from 
the nematic superfluid to an isotropic normal state is pos- 
sible. These predictions are straightforward to test exper- 
imentally in ultracold atomic g vortex-unbinding 
can be seen via interferometry |42] and orientational or- 
der via time-of- flight imaging [43]. The nature of the 
interplay between dislocations, disclinations, and exotic 
bound states can be generalized from superfluidity to 
e.g. spin-density-wave order as in the condensation of 
magnons with spin-orbit coupling |44j . Finally, we have 
argued that the spatial extent of the modes bound to 
dislocations in weak topological insulators or supercon- 
ductors is governed by the disclination size (i.e., the in- 
verse orientational stiffness), which is typically smaller 
than the dislocation size (i.e., the inverse shear stiffness). 
The extensions of this classification of disclinations to 
other classes of states, such as nematic Chern insula- 
tors [J51 SB], and to quantum melting [171 EE] will be 
considered in future work. 
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